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Introduction 

Hydrogen  is  an  important  technological  component  for  silicon  surface  passivation  during 
the  fabrication  of  metal  oxide  semiconductor  (MOS)  devices.  In  a recent  development 
[1]  deuterium  has  been  substituted  for  hydrogen  in  the  process  of  passivation.  This 
resulted  in  the  significant  decrease  of  hot  carrier  degradation.  Being  isotopes,  deuterium 
and  hydrogen  have  vitually  identical  electronic  properties.  Therefore  the  significant 
difference  in  adatom  behavior  as  observed  in  [1]  must  be  related  to  the  surface  lattice 
dynamics. 

The  significant  difference  in  the  isotope  atoms  behavior  occurs  not  only  under  ex- 
treme nonequilibrium  conditions  but  even  in  the  equilibrium  properties  of  isotope  atoms 
on  the  surface. 

This  paper  deals  with  the  calculations  of  the  concentrations  of  adatoms  H (or  D)  on 
the  surface  of  Si  which  is  in  thermodynamic  equilibrium  with  hydrogen  or  deuterium 
gas. 

'Hie  description  of  the  equilibrium  between  the  Si  surface  and  the  gas  of  molecules 
requires  the  knowledge  of  both  the  chemical  potential  of  H(D)  on  the  surface  and  the 
chemical  potential  of  H(D)  in  the  gas.  Along  with  the  electronic  contribution,  the  chem- 
ical potential  of  the  adatom  has  the  vibration  contribution  which  can  be  found  from  the 
lattice  dynamics  of  the  impurity  atoms  on  the  surface.  'Hie  chemical  potential  of  H(D) 
atom  in  the  gas  is  known  from  the  thermodynamics  of  two  atomic  gases.  The  equality 
of  the  chemical  potentials  of  H(D)  atoms  on  the  surface  and  in  the  gas  in  the  ther- 
modynamic equilibrium  enable  us  to  find  the  surface  equilibrium  concentrations  of  the 
adatoms.  The  difference  of  H and  D surface  concentrations  results  from  the  difference 
of  H and  D surface  localized  vibrations.  It  is  shown  that  the  surface  concentration  of 
D-adatoms  is  by  10  to  15  times  higher  then  the  concentration  of  H-adatoms. 

1 Dynamics  of  adatoms  on  the  surface  of  Si 

In  order  to  find  the  surface  localized  vibration  frequencies,  we  consider  the  general 
form  of  the  equation  of  the  motion  for  the  semiinfinite  crystal  with  some  adatoms 


(Z  — D)u  = 0 (1) 

where  Z represents  vibrations  of  the  perfect  crystal  with  the  surface  and  D is  the 
perturbation  matrix  of  the  adatoms,  u is  the  column  of  the  displacement  of  all  the 
atoms. 

For  simplicity  the  nonreconstructed  surface  (100)  of  Si  with  the  symmetry  of  the 
simple  square  is  considered.  We  assume  that  the  adatom  of  the  mass  M(H)  is  bound 
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in  the  “on-site”  configuration  to  the  surface  atom  of  Si  with  the  mass  M by  the  force 
constant  7 /.  Restricting  ourself  by  the  nearest  neighbour  approximation  for  the  descrip- 
tion of  the  interaction  of  the  adatoms  with  the  crystal,  one  can  get  for  the  case  of  the 
low  concentration  of  the  adatoms  the  following  characteristic  equation  for  finding  the 
adatom  localized  vibration  frequency 

1 - -4^G(<J)  = 0,  (2) 

LOz  — LO  Q 

where  cUq  = 7 ' /M^H\  G(ui2)  is  the  diagonal  element  of  the  Green’s  function  matrix 
G = l-\ 

The  Green’s  function  G(co2)  can  be  given  in  terms  of  the  frequency  distribution 
function  of  the  perfect  crystal  go(to2)  [2] 


, 1 r“i 

ReG(“)  = s/o 


go(Q2) 

tu2  - Q2 


dCl2. 


(3) 


Here  go  (to2)  is  normalized  to  the  unity.  The  quantity  u>L  is  the  maximum  frequency  of 
the  crystal  with  the  free  surface. 

Since  the  localized  vibration  frequency  of  the  light  adatom  tuioc  satisfying  to  Eq.  (2) 
is  significantly  higher  than  the  maximum  frequency  of  the  Si  phonon  spectrum,  the 
main  contribution  in  the  integrant  at  to  = iuioc  is  given  by  the  high  frequency  optical 
vibrations.  Because  the  optical  branches  have  the  small  dispersion,  the  important  optical 
part  of  the  distribution  function  could  be  approximated  by  a 5-function 

go(to2)  « goP‘(y2)  = S(uj2  - u>l).  (4) 


From  Eqs.  (2)-(4)  the  algebraic  biquadratic  equation  for  the  determination  of  the 
surface  localized  frequency  tuioc  is  obtained: 


1 - 


to 


1 (u)2  + ) 


(tO2  — £Uq)<U2 


= 0. 


(5) 


Using  the  numerical  value  of  the  force  constant  7/  from  [3]  7/  = 8.8  x 104din  • cm  1 , 
we  calculated  from  Eq.  (5)  the  following  localized  frequencies  of  H (or  D)  stretching 
localized  modes:  cu®  = 2.34  x 1014  s_1  = 1.70  x 1014  s-1.  We  note  that  cu®  > 

uC . These  frequencies  differ  from  the  experimental  data  cu®  = 3.96  x 1014s~~‘, 
<*4^  = 2.83  x 1014s-1  [4],  The  discrepancy  is  apparently  the  result  of  our  simple 
approximation. 

It  is  shown  in  [5]  that  the  account  of  the  next  nearest  neighbour  approximation 
and  of  the  noncentral  forces  would  result  in  bending  modes  obtained  along  with  above 
mentioned  stretching  modes.  The  experimental  value  of  the  double  degenerate  at  T — 
point  bending  mode  of  H(D)  on  the  surface  of  Si  is  cu®d  = 1 . 17  x 1014  s^1  [4],  Because 
the  bending  mode  frequency  is  lower  then  the  stretching  mode  frequency,  stretching 
modes  are  the  most  important  for  the  calculation  of  the  partition  function  of  the  system. 
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2 Chemical  potential  of  H(D)  adatom  on  the  surface 


In  our  case  of  the  low  surface  concentration  of  adatoms  they  constitute  a two-dimensi- 
onal dilute  solution  on  the  surface  of  Si.  The  vibrational  part  of  the  partition  function 
of  the  solution  can  be  given  in  terms  of  the  frequency  distribution  function  g(ui2)  of 
the  crystal  with  adatoms  on  the  surface 


Zcrvst  — 


Nmifl 


WadatomKAWf  - ^adatom)! 


exp 


(6) 


Here  7VSUIf  is  the  total  number  of  the  dangling  bonds  on  the  surface  and  Adatom  is  the 
number  of  the  adatoms.  The  combinatorial  factor  represents  the  possibility  of  inter- 
changing within  both  the  occupied  and  empty  surface  sites.  The  frequency  cumax  is  the 
maximum  frequency  of  the  vibrations  for  the  crystal  with  adatoms  on  the  surface. 

With  high  accuracy,  the  total  frequency  distribution  function  may  be  given  by 

g(oj2)  = 3A/sigo(w2)  + Nadatom36(u>2  - cu,20C).  (7) 


Substituting  (7)  in  (6)  and  using  the  equation  for  Helmgolz  free  energy  F = —Tin Z 
and  the  equation  for  the  chemical  potential  of  the  adatom  ^adatom  = {dF  / dNadat0m)TV, 
we  get  the  chemical  potential  for  hydrogen  adatom 


Ms”rf  = nn«idatom  + 3Tln 


2 sinh 


hw. 


(H) " 

loc 


2 T 


and  for  deuterium  adatom 


/'2  = -i  rn" 


2si„h?^ 


2 T 


+ £h> 


+ £d- 


(8) 

(9) 


Because  binding  energies  of  isotopes  are  equal,  £h  = £d  and  the  difference  of 
chemical  potentials  comes  from  the  difference  of  localized  frequencies  ^ oj'^J 
only. 


3 Chemical  potential  of  H(D)-atoms  in  the  gas  phase 


The  chemical  potentials  of  H and  D atoms  in  the  gas  is  obtained  from  the  conditions  of 
the  chemical  equilibrium  for  the  reaction  of  ^(Dj)  dissociation  and  from  the  chemical 
potential  of  molecule  consisting  of  two  identical  atoms  [6]. 


(H)  = 

gas 


1 


[rin/^2  -CpTlnT  - Ch 2T  +e0] 


(10) 


and 

l42  = \ \T  ^ Pd2  - cP  T In  T - Cm  T + e0]  ( 1 1 ) 

Here  /,H2('f>D2)  is  the  gas  pressure,  cp  is  the  specific  capacity  at  the  constant  pressure, 
Ch2  = frl[(7H2/ft5)(M^2/27^)3/,2]  and  Cd2  = hi[(/D2/7t5)/(A7D2/27r)3/2]  are  the  chemical 
constants  of  the  hydrogen  gas  or  the  deuterium  gas,  respectively  (/h2  and  7D2  being  the 
momentum  of  inertia  of  the  molecules),  eo  is  the  binding  energy  of  the  molecule  in  the 
gas. 
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4 Concentration  of  H(D)  atoms 

The  substitution  of  chemical  potentials  (8)  and  (10)  for  hydrogen  or  (9)  and  (11)  for 
deuterium  in  the  condition  of  thermodynamic  equilibrium  /x®  = p^s)  and  /4urf  = 

CD’) 

/igas  , respectively,  allows  to  find  the  concentration  of  hydrogen  atoms  and  the  concen- 
tration of  deuterium  atoms  on  the  Si  crystal  surface  in  the  thermodynamic  equilibrium 


n 


(H) 

adatom 


n 


(D) 

adatom 


pV2Tc„/2 


pV2Tc?/2 


( M\\2' 

3/2' 

1/2 

[ /Tl 

exp 

j£q  — £h 

V 27r , 

) 

exp 

T 

T 

h)7  f 
h 5 I 27T 


3/2' 


1/2 


exp 


(D) 


exp 


Je0  ~ £D 


(12) 

(13) 


Since  the  localized  frequencies  tul^)  and  satisfy  the  condition  T < 
and  /h2//d2  = Mh2/Mv2,  one  can  get  the  following  equation  for  the  ratio  of  the  con- 
centrations 


n 


(H) 

surf  


n 


(D) 

surf 


/MhA5/4 

\Md2J 


exp 


3^[4!')  - 

2 T 


(14) 


It  is  seen  from  Eq.  (14)  that  the  concentration  of  deuterium  exceeds  significantly 
the  concentration  of  hydrogen  due  to  the  difference  of  both  the  isotope  masses  and  the 
localized  frequencies. 

Using  experimental  values  for  localized  frequencies  of  H and  D we  have  found 
that  the  ratio  of  concentrations  is  about  15  at  the  typical  technological  temperature 
T = 700  K.. 
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